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ABSTRACT 

In recent years, detailed observations and accurate numerical simulations have provided support to the idea 
that mergers of compact binaries containing either two neutron stars (NSs) or an NS and a black hole (BH) 
may constitute the central engine of short gamma-ray bursts (SGRBs). The merger of such compact binaries 
is expected to lead to the production of a spinning BH surrounded by an accreting torus. Several mechanisms 
can extract energy from this system and power the SGRBs. Here we connect observations and numerical 
simulations of compact binary mergers, and use the current sample of SGRBs with measured energies to 
constrain the mass of their powering tori. By comparing the masses of the tori with the results of fully general- 
relativistic simulations, we are able to infer the properties of the binary progenitors which yield SGRBs. By 
assuming a constant efficiency in converting torus mass into jet energy, ejct = 10%, we find that most of the tori 
have masses smaller than O.Olil/0, favoring "high-mass" binary NSs mergers, i.e., binaries with total masses 
> 1.5 the maximum mass of an isolated NS. This has important consequences for the gravitational-wave 
signals that may be detected in association with SGRBs, since "high-mass" systems do not form a long-lived 
hypermassive NS after the merger. While NS-BH systems cannot be excluded to be the engine of at least some 
of the SGRBs, the BH would need to have an initial spin of ~ 0.9, or higher. 

Keywords: accretion, accretion disks — gamma-ray burst: general — gravitational waves — methods: numer- 
ical — stars: neutron 



1. EvlTRODUCTION 

Binary neutron star (BNS) and neutron star-black hole (NS- 
BH) binaries are the leading candidate s for the central engine 
of short gamma-ray bursts (SGRB s; iBlinnikov et al.l [T984t 
iPaczyn ski 1986*: lEichler et al.l ll989^. They are also one of the 
most powerful sources of gravitational waves (GWs), and ad- 
vanced interferometric detectors are expected to observe these 
sources at rates of 0.4 — 400 and ^ 0.2 — 300 events p er 
year for BNS and NS-BH, respectively jAbadie et al.ll201(il) . 

Fully general-relativistic simulations have shown how 
such mergers can lea d to the formation of accretion disks 
around spinning BHs (iBaiotti et al.l 

o| 



iKiuchi etaPIIOOgt IFaber & Rasi 



2008 



lEtienne et alJl2009l: 



2012 ). Moreover, when 



magnetic fields are present, they can provide one of the mech- 
anisms neces sary to extrac t energy, and power collima ted rel- 
ativistic jets ( iRezzolla et a l. 2011; Etienne et al.ll2012l) . 

So far, properties of the progenitors of SGRBs have 
been inferred by studying their re dshift distribution , 
close environmen t, host galaxy types (iBloom et al.l 120061 : 
IZhang et alJ l2009t lBergeili201li and by comparing those 
observat ions with predictions from population synthesis 
models (Perna & Belczvnski' 120021; iBelczvnski etall 120061; 
lO'Shaughnessy et al. 2008). In this Letter we make a con- 
nection between theory and observations, which allows us 
to directly probe the SGRB progenitors. In particular, we 
consider a complete (to date) sample of SGRBs with mea- 
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sured redshifts, and Unk the properties of their observed emis- 
sions to the masses of the tori responsible for their genera- 
tio n. By comparing thes e tori w ith the theoretic al predictions 
of IRezzolla etaP (l20Toh and of lFoucartI JMl) . we are able 
to infer the properties of the compact binaries that may have 
generated such bursts. 

In Section |2] we provide details on the sample of SGRBs 
considered in this Letter In Section[3] we use the theoretical 
results to compute the masses of the tori that have generated 
such bursts, and link them to their progenitors. In Section |4] 
we show how GWs may be used to further constrain the pro- 
genitors, and in Section|5]we summarize our main results. 

2. GRB SAMPLE DATA 

We selected our sample of SGRBs based on three crite- 
ria: duration, hardness ratio, and spectral lags. Swift SGRBs 
with known redshift are listed in Table [T] SGRBs with a 
temporally extended emission (EE) were also considered. In 
the latter case, the quoted energetics include the contribu- 
tion of the short-hard spike, and of the EE. Since the two 
emi ssion episodes typ ically have a comparable energy bud- 
get jNorris et alj|201 lb . the presence of EE affects our calcu- 
lations by a factor ss 2. 

The burst energetics, iS^- iso (Column 3 in Table [l), were 
calculated by using the prornpt ern issi on spectral param eters 
(mainly from lSakamoto et aLlEoTll and lGoldst ein et al."2012|) 
and sh ifted to a common rest-frame energy band (.Bloom et aP 
I2OOI . When possible, we used measurements of the broad- 
band GRB spectrum (e.g., by the FennUGBM) and calculated 
i?-y_iso in the comoving 10 keV-10 MeV energy range. In most 
cases, only SwiftfBAI observations are available, and we re- 
port the burst energetics in the narrower 15-150 keV rest- 
frame band, thus unavoidably underestimating th e bolomet- 
ric en ergy release. For a typi cal Band spectrum ( iBand et al] 
11993b . peaking at « 500 keV jNavaet al.ll2011l) . we estimate 
an average /c-correction factor of 6. We therefore do not 
expect that the uncertainty in the GRB spectral shape of Swift 
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Table 1 

SGRB Sample 



Table 2 

BNS Simulations and Toms Masses 



GRB Name 


z 






1 (erg) 


AA' (keV) 




rue 


■AMq) 


050509B 


0.225 


9.1 


X 


10*^ 


15 - 


150 


1.0 


X 


10~ 


-b 


050709(EE) 


0.161 


3.4 


X 


10« 


10 — 


10^ 


3.8 


X 


10~ 


-4 


050724(EE) 


0.257 


1.9 


X 


1050 


15 — 


150 


2.1 


X 


10" 


■3 


051221A 


0.546 


2.9 


X 


1051 


10 — 


10"* 


3.3 


X 


10" 


.2 


061006(EE) 


0.438 


2.1 


X 


1051 


10 — 


10"* 


2.4 


X 


10" 


-2 


070429B 


0.902 


2.1 


X 




15 — 


150 


2.3 


X 


10" 


-3 


070714B(EE) 


0.923 


1.6 


X 


1052 


10 — 


10"* 


1.8 


X 


10" 


-1 


071227(EE) 


0.381 


1.2 


X 


1051 


10 — 


10"* 


1.4 


X 


10" 


-2 


080905A 


0.122 


4.5 


X 


10*» 


10 — 


lO" 


5.1 


X 


10" 


-4 


090510 


0.903 


4.7 


X 


1052 


10 — 


10* 


5.2 


X 


10" 


-1 


1001 17A 


0.920 


1.4 


X 


1051 


10 — 


lO'i 


1.6 


X 


10" 


-2 


111117A 


1.3 


5.3 


X 


1051 


10 ~ 


10"* 


6.0 


X 


10" 


-2 


051210 


1.3 


4.0 


X 


1050 


15 - 


150 


4.5 


X 


10" 


-3 


060801 


1.130 


1.9 


X 


1050 


15 - 


150 


2.1 


X 


10" 


3 


061210(EE) 


0.410 


5.6 


X 


1050 


15 - 


150 


6.2 


X 


10" 


■3 


070724A 


0.457 


2.3 


X 


10« 


15 - 


150 


2.5 


X 


10" 


-4 


070729 


0.8 


1.6 


X 


1050 


15 - 


150 


1.8 


X 


10" 


■3 


080123(EE) 


0.495 


5.7 


X 


1050 


15 - 


150 


6.3 


X 


10" 


3 


101219A 


0.718 


7.4 


X 


1051 


10 - 


10" 


8.2 


X 


10" 


-2 


060502B 


0.287 


9.8 


X 




15 - 


150 


1.1 


X 


10" 


-4 


061217 


0.827 


6.8 


X 


10« 


15 - 


150 


7.6 


X 


10" 


-4 


061201 


0.111 


9.4 


X 


10« 


15 - 


150 


1.1 


X 


10" 


-4 


070809 


0.473 


7.9 


X 


10« 


15 - 


150 


8.8 


X 


10" 


-4 


090515 


0.403 


1.0 


X 


10« 


15 - 


150 


1.2 


X 


10" 


-4 



Note. — The different columns refer respectively to the GRB name, the 
redshift z derived from the GRB host, the isotropic equivalent gamma-ray 
energy -E-,,iso. measured in the rest-frame energy band AS, and the mass 
of the torus A/torus (see Equation (TJ). The different blocks ref er to the 
uncertainty in the SGRB/host galaxy association IBloom et alJ2002l) . The top 
one includes SGRBs with a precise identification of a host galaxy; those in 
the middle have a less certain association with their host; those in the bottom 
are significantly offset from the associated host galaxy, and are affected by a 
larger uncertainty. 



bursts may have a major impact on the resuhs. 

Table [T] shows that SGRBs display a wide range of en- 
ergies, from 10*^ erg to 10^^ erg, with a median value of 
2 X 10^° erg. The quoted values refer to the isotropic equiv- 
alent gamma-ray energy, while the true energy scale also de- 
pends on the outflow beaming factor /h = 1 — cos(6'jct), be- 
ing 6'jot the jet opening angle. The degree of collimation of 
SGRBs is still a poorly constrained quantit y, inferred values 
range from k. 0.001 to « 0.1 (Bu rrows et al.l 120061: 
iNicuesa Guelbenzu et alj|2012h . but in most cases only weak 
lower bounds can be placed. Here, we use the isotropic ener- 
gies listed in Table [T] to set an upper limit to the burst-energy 
release. 

Assuming that all the SGRBs in our sample were produced 
by accretion tori around spinning BHs, we now correlate the 
values for the isotropic energy listed in Table [T] with the mass 
of such tori. In particular the torus mass is determined as 



eMu 



(1) 



where e is the efficiency in converting the mass of the torus 



into the isotropic gamma-ray emission E, 



7,iso- 



Here, 



e is given by the product of two efficiencies: one to convert 
mass of the torus into jet energy, ejot, and the other to convert 
the latter into gamma-rays, e^. 

Fully general-relativistic simulations of BNS mergers 
have shown the formation of thi n and highly mag- 
netize d tori around spinning BHs jRezzolla et al.l 120101 
1201 lb . Here we make the important assumption that 
SGRBs are powered via magnetic fields (.Blandford & Znaie3 



Model 


A^BNS 


q 


-^torus 




A/SNS/A'^max 




(A^0) 




(Mq) 


(Mq) 




1 .46-45-IF 


3.24 


1.00 


0.1374 


2.20 


1.47 


1 . 62-45-IF 


3.61 


1.00 


0.1101 


2.20 


1.64 


M3 . 6ql . 00 


3.90 


1.00 


0.0012 


2.20 


1.77 


M3 . 7qO . 94 


4.03 


0.94 


0.0121 


2.20 


1.83 


M3 . 4qO . 91 


3.76 


0.92 


0.1202 


2.20 


1.71 


M3 . 4qO . 80 


3.72 


0.81 


0.2524 


2.20 


1.69 


M3 . 5qO . 75 


3.80 


0.77 


0.1939 


2.20 


1.73 


M3 . 4qO .70 


3.71 


0.72 


0.2558 


2.20 


1.69 


APR145145 


2.87 


1.00 


0.000549 


2.18 


1.32 


APR1515 


2.97 


1.00 


0.000134 


2.18 


1.36 


APR1316 


2.87 


0.81 


0.0275 


2.18 


1.32 


APR135165 


2.97 


0.82 


0.00707 


2.18 


1.36 


APR4-2 8 


2.77 


1.00 


0.003 


2.21 


1.25 


SLy-27 


2.67 


1.00 


0.02 


2.05 


1.30 


H3-27 


2.68 


1.00 


0.05 


1.79 


1.50 


H3-29 


2.87 


1.00 


0.01 


1.79 


1.61 


H4-27 


2.68 


1.00 


0.18 


2.03 


1.32 


H4-29 


2.87 


1.00 


0.02 


2.03 


1.41 


H4-30 


2.97 


1.00 


0.01 


2.03 


1.46 


ALF2-27 


2.67 


1.00 


0.16 


2.09 


1.28 


ALF2-29 


2.87 


1.00 


0.02 


2.09 


1.38 


ALF2-30 


2.97 


1.00 


0.003 


2.09 


1.42 


PS-27 


2.68 


1.00 


0.04 


1.76 


1.53 


PS-29 


2.88 


1.00 


0.02 


1.76 


1.64 


PS-30 


2.97 


1.00 


0.01 


1.76 


1.69 



Note. — The different columns represent respectively the name of the 
model, the gravitational mass of the binary, A/bns> the mass ratio of the grav- 
itational masses of the two NSs, q, the baryonic mass of the torus. A/torus, 
the maximum gravitational mass of an isolated NS for the equation of state 
(EOS) used in that simulation. A/max, and the ratio between the mass of the 
binary and A/^ax, A/sNs/^max- Th e different blocks of the table refer, 
from top to bottom, to the simulations by'Baiotti et al.' 1*20081. "Rezzo lla et al.l 
t2010), Kiuchi et al. (2009), and Hotokezaka et al. (201 1). Note that the sim- 
ulations reported in Bai otti et alj ( ,200 8i) and iRezzolla et al.. t2010i) used an 
ideal-fluid EOS and hence they can be scaled to different masses. Here we 
have chosen the values for an ideal- fluid EOS so that Af max = 2.20, in 
agreement with current observations jPemorest et aU2010l) . 



[19771 iBlandford & Pavria [T9821: iTchekhovskovetaD lIoTl) 
and ignore the effects of vis cosity and neutrino cool- 
ing dChen & Beloboro dov 2007')FI General-relativistic mag- 
netohydrodynamic (GRMHD) simulations of accretion disks 
showed that the efficiency in converting torus mass and 
BH spin into jet energy (i.e., ejet) varies between few 
per cent up to more than 100% for maximally spin- 
ning BHs ('De Villiers et al."2005': 'Tche khovskov et al.l l201 It 
iMcKinney et al. 2012: Fragile et al. 20121). The efficiency de- 
pends sensitively on the BH spin, the disk thickness, and the 
magnetic flux. Accounting for all of these effects is currently 
not possible, and thus we made the simplifying assumption 
of a constant efficiency for all BNS and NS-BH mergers, 
ejct = 10%. While an obvious approximation, our main re- 
sults do not change sensitively if ejot is taken to be larger than 
-0.1%. 

After the jet is emitted, a fraction of its energy is converted 
into gamma rays. The conversion efficiency f or a sample 
of lon g and short Swift GRBs was computed bv lZhang et alj 
(l2007l) bv comparing the gamma-ray fluence with the bright- 
ness of the X-ray afterglow at early and late times. They find 
that, while the efficiency in long GRBs varies strongly from 
burst to burst, ranging from a fraction of a per cent to almost 
100%, in SGRBs the range is narrower, varying between 30% 

^ For a discussion of neutiino -powered SGRB s and the relation betwee n 
simulations and observations, see lLee et al.U200^ . IOechslin & Jan"k3 120061) . 
and Fan &Wei( 20113. 
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Figure 1. Left panel: MtoTue/MQ as a function of the mass ratio q and of the ratio between the gravitational mass of the binary and the maximum mass for 
an isolated NS (Mbns /Afmax). Aftorus has been computed using Equation j2). The dotted lines are the isocontours corresponding to the Mtorus values in 
Table[T] Right panel: plot of Mtorus as a function of A/eNs/A^max for all the equal-mass (q = 1) simulations reported in Table|2] The horizontal bars give the 
percentage of the SGRBs in Table[T]that are generated by tori with that range of masses (assuming a total efficiency e of 5% as in Table[T). Since the mass of the 
torus increases for g < 1 . each point should be considered as a lower limit on the mass that can be obtained for that EOS and mass of the binary. The different 
points refer to A/torus computed from simulations of BNS mergers using different EOSs (see Table|2j- 



and 60%, with an average of 49%. We hence assume a fiducial 
value of = 50%, so that the total efficiency in Equation ([T]i 
becomes e = 5%. The last column of Table [T] shows the cor- 
responding torus masses. 

3. TORUS MASSES 

In th e following, we lin^A/torus to the theoretic al predic- 
tions of lRezzoUa et all (120101) and of lFoucarj (120121) . who de- 
rived analytic fits from the results of fully general-relativistic 
sim ulations of BNS and NS-BH mergers, respectively (see 
also lPannaraleetal.ll2011h . 

3.1. Binary Neutron Star Mergers 

iRezzolla et al.l (1201 Oh derived a phenomenological expres- 
sion to compute the masses of the tori formed by BNS merg- 
ers. Here we have revised that fit and expressed it as a function 
of two dimensionless quantities: the gravitational mass ratio 
q < 1 and the ratio between the gravitational mass of the bi- 
nary and the maximum gravitational mass for an isolated NS 
(MBNs/Mmax). We derived 



[Ci(l -q)+ C2][C3(1 +q)- A/BNsMnax] • (2) 



The coefficients Cl = 2.974±3.366, C2 = 0.11851±0.07192, 
and C3 = 1.1193 ± 0.1579 were determined by fitting Equa- 
tion (|2]i t o the results of the full y general-relativistic s imu- 
lations of iBaiotti etlTI ( 120081) and IRezzolla et all (1201 Oh . but 
rescaled to allow for a value of M,„ax = 2.20 Mq to be more 
consistent with current observations of NS masses (see also 
TableEJ. 

The left panel of Figure [T] shows Mtorus computed using 
Equation (Q as a function of q and Mbns /Mmax, while each 
of the red dotted lines represents the isocontour relative to an 
observed GRB in Table [T] when assuming our fiducial value 
e = 5%. The right panel shows the distribution of torus 
masses obtained from observations (horizontal bars, see Ta- 
ble [TJ together with the mass of the torus computed from nu- 
merical simulations of equal-mass BNSs {q = 1, see Table[2]i. 
As one can easily see, two thirds of the SGRBs of our sam- 
ple appear to be generated by tori with masses smaller than 
~ lO^^M©. Moreover, since A/torus increascs for 5 < 1, 



each point should be considered as a lower limit on the mass 
that can be obtained for that equation of state (EOS) and 
BNS mass. This means that, while the energetics of most 
SGRBs can be explained by current numerical simulations, 
some of the less energetic SGRBs should result from BNS 
mergers with masses larger than the ones simulated so far 
(since A/torus dccrcascs with increasing .A/bns)- 

It is evident from the left panel of Figure[T]that, for an ideal- 
fluid EOS, almost all of the SGRBs would be generated by 
BNSs with A/BNs/-^max > 1-8 and hence they would be 
"high-mass" systems. This means that the mass of the sys- 
tem would be too high to lead to the formation of a long-lived 
hypermassive NS (HMNS) and that the merger would pro- 
duce a prompt collapse to BH. This is also true for the mod- 
els with realistic EOSs shown in the right panel of Figure [T] 
For example, the two circles refer to simulations of equal- 
mass binaries using an APR EOS (models APR14 514 5 and 
APRl 5 1 5; see lKiuchi et al.ll2009h and they produce tori with 
masses in the rang e of ~ 33% of all SGRBs in our sam- 
ple. As reported in iKiuchi et al.l (120091) . in both cases, col- 
lapse to BH occurs 1 ms after merger. We recall that the 
threshold A/bns /M^ax below which a long-lived HMNS is 
formed is strongly dependent on the EOS. All the simula- 
tions that produce tori with masses < 0.1 Mq in the right 
panel of Figure [T| produce an HMNS that collapses on a 
timescale of few ms jKiuchi et al.ll2009t IRezzolla et al.ll20Tol: 
iHotokezaka et al.l20I ll) PI As we discuss in SectionjH this has 
a fundamental impact on the GW signal we may expect from 
SGRBs. 

3.2. NS-BH Mergers 

iFoucarj (120 12h derived the following fit for the mass of the 
torus produced by an NS-BH merger: 



(M^ 



NS 



\M^s 



^(3/(7)1/3(1 - 2Cns)Mns - fiRiscoCns 

(3) 



If e was one order of magnitude smaller, Mtorus would be 10 times 
larger, but 67% of the SGRBs would still have A/torus < 0.1 A/q and hence 
be generated by "high-mass" systems. 
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Figure 2. Similai' to the left panel of Figure [T] but for the NS-BH case. A/torus is shown as a function of the NS compactness Cns and mass ratio 1/g = 
orus is computed using Equation {sj. Each panel assumes a different value for the dimensionless spin of the BH, x- In all cases, we assume 
A^NS = 1.4 Af0. The dotted lines are the isocontours corresponding to the A/torus values in Table[T] 



where a = 0.288 ± 0.011, [3 ^ 0.148 ± 0.007, M^s is the 
gravitational mass of the NS, its baryonic mass, Cns the 
NS compactness, l/q = A/bh/A/ns > 1 the ratio between 
the BH and NS mass es, and J^isco th e radius of the innermost 
stable circular orbit dFoucarli 120121) . We note that in order 
to compute Mtorus we need to know the ratio Mj^s/M^c,, 
for which there is no analytic expression available. We make 
here the reasonable assumption that the baryonic mass is 10% 
larger than the gravitational massQ This assumption may lead 
to a few percent error on the mass of the torus, which is suffi- 
ciently small to not affect the results of this Letter 

The four panels in Figure |2] show A/toms, computed 
using Equation (|3]l, as a function of the NS compactness 
Cns and mass ratio l/q. Each panel assumes a different 
value for the dimensionless spin of the BH, x- It is evident 
from this figure that not even the most rapidly spinning 
BH (x = 0.9) can explain the most energetic burst in our 
sample (GRB090510 with A/torus ^ O.SMq). Moreover, if 
we acc ount for the results of populations synthesis calcu- 
lations (iBelczvnski et al.l 120081) . which predict most of the 
NS-BH binaries to have mass ratios l/q between ^ 7 and 
~ 10, w hile the NS compac tness is expected to be larger than 
- 0.16 dSteiner et al.ll2012l) . then most of the SGRBs in our 
sample can be only explained if the binary has a BH with 



For the NSs reported in Table|2] the baryonic mass A'/£g is ~ 8% larg er 
than the gravitational mass A/^jg for the ideal-fluid EOS and ~ 11% larger 
for the APR EOS. 



an initial spin of ^ 0.9 or larger From current observations 
of SGRBs, it is then clear that, while current simulations of 
BNS mergers may easily produce tori in the range required to 
explain all the current observations, NS-BH mergers cannot 
be used to explain the most energetic bursts. 



4. CONSTRAINTS USING FUTURE GW OBSERVATIONS 

As shown in lBaiotti et"an (120081) and lRezzoUa et alj (l2010t) . 

the GW signal is strongly affected by the mass of the sys- 
tem and how close this is to the maximum mass for each 
particular EOS. BNSs with masses close to the maximum 
mass exhibit a prompt collapse to BH after the merger, 
while lower-mass systems produce an H MNS which can sur- 
vive fro m few ms up to hundreds of ms ("Ba iotti et alj|2008l : 
iRezzoll a et al. 2010; Giacomazzo et al. 201 ijl GW signals 
from "high-mass" systems are simply composed of the in- 
spiral, merger, and BH ring-down phases. Lower-mass sys- 
tems, instead, display a more complex GW signal with a 
rich spectrum due to the emission of GWs from the HMNS 
formed after the merger Such emission is i mportant since it 
can h elp infer the properties of the NS EOS (Ba uswein et alJ 
l20T2h . However, since the GW signal emitted by the 
HMNS i s in a range of frequencies between ~ 2kHz and 
~ 4kHz ( iBaiotti et alJl2008HBauswein et alj|2012h , it may be 
difficult for advanced LIGO/Virgo to detect it, and a third gen- 
eration of detectors, such as the Einstein Telescope, would be 
required dAndersson et al...201 1,) . On the other hand, the for- 



Progenitors of Short Gamma-Ray Bursts 



GRB 080905A 



GRB 050724 



GRB 07071 4B 



A* 



= 10"'' . ♦ ■ T 



1.4 



1.6 

/-^' 



1.2 1.4 1.6 



1.4 1.6 

Mam /."„..,« 



Figure 3. Different panels plot A/torus as a function of A/eNs/A/max for all the simulations reported in Table|2] and compare them with three SGRBs taken 
from TableH GRB080905A, GRB()50724, and GRB070714B. In each panel, a horizontal dashed line represents the value of A/torus reported in Tablenwhile 
the shaded region represents the range of A/torus assuming a total efficiency between e = 1% and e = 10%. Symbols for the various EOS from Table|2]are the 
same as in Figure[T] but here also the q 1 simulations have been included. 



mation of an HMNS after the merger can also be inferred by 
measuring the delay time between the BNS merger (indicated 
by the GW signal) and the time of the emission of the SGRB 
(which we may assume coincident with BH formation). A 
delay time of ^ 100 ms or larger would clearly indicate the 
formation of an HMNS. 

As discussed in Section lTTl we find that only the most en- 
ergetic SGRBs can be compatible with a low-mass binary and 
hence the formation of an HMNS. For the greatest majority of 
SGRBs, a high-mass system is the most likely scenario, and 
hence we expect SGRBs to be observed simultaneously with 
GWs which would lack the high-frequency emission typical 
of the HMNS. Although the GW signal from a prompt col- 
lapse is not as rich as that from an HMNS, the simultaneous 
detection of an SGRB with the associated GW may help con- 
siderably in constraining the NS EOS. 

This is illustrated in Figure [3] with a selection of three 
SGRBs from our sample (with low, medium, and high en- 
ergetics). In each panel, the horizontal dashed line represents 
the value of A/toms reported in Table [T] while the shaded re- 
gion represents the range of il/torus assuming a total efficiency 
e between 1% and 10%. The various points represent A/torus 
computed from simulations of BNS mergers using different 
EOSs (see Table |2]i. If an SGRB was detected together with 
a GW signal, we could use the energetic of the burst to de- 
termine the torus mass (horizontal bars in FigureO, while the 
GW signal could be used to infer the mass of the BNSs (which 
would give a vertical bar in those panels). The combination 
of these two pieces of information would restrict the allowed 
EOS parameter space. 

In the case of NS-BH binaries, the mass ratio q and the spin 
of the BH can in principle be measured via GW observations; 
then a simultaneous detection of a GW and an SGRB would 
allow to set an independent constraint on the NS compactness 
Cns and hence infer the NS EOS. 

5. SUMMARY 

We have performed a novel analysis of the energetics of 
SGRBs in connection with the properties of the compact bi- 
nary systems that may have generated them. We have shown 
that most of the SGRBs could be produced by magnetized 
tori with masses lower than ^ O.OIMq. Combining this in- 
formation with the results of numerical simulations of NS- 
NS and NS-BH mergers, we have concluded that most of the 
SGRBs are consistent with the merger of "high-mass" BNS 
systems, i.e., with MeNs/Afmax > 1-5. While NS-BH sys- 
tems cannot be completely excluded, the BH would need to 
have an initial spin of ^ 0.9, or higher. Moreover, the most 



energetic SGRBs, such as GRB090510, could not be pro- 
duced by the merger of an NS with a BH[[3 We note that 
while our results are affected by some uncertainty in the exact 
value of the efficiency e, our conclusions are robust as long as 
ejot ^ 0.1% (i.e., e ^ 5 X 10"''), which is much lower than 
what was observed in GRMHD simulations of jets from accre- 
tion disks (De Villiers gt_ al. 2005; Tchekhovskoy et al. 2011]; 
iMcKinnev et al.ll2012tlFfagile et al.ll2012l) . 

GW signals from SGRBs would help validate our results. In 
particular, in the case of BNSs, since we find that SGRBs are 
most likely generated by "high-mass" BNSs, the GW signal 
would lack the features that are associated with the formation 
of an HMNS, since "high-mass" BNSs produce a prompt col- 
lapse to BH a few ms after merger. A simultaneous detection 
of GWs with SGRBs would help constrain the EOS of NS 
matter. 
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